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ABSTRACT 

We present a high-resolution (R^ 50 000) atlas of a uranium-neon (U/Ne) hollow- 
cathode spectrum in the H-band (1454 nm to 1638 nm) for the calibration of near- 
infrared spectrographs. We obtained this U/Ne spectrum simultaneously with a laser- 
frequency comb spectrum, which we used to provide a first-order calibration to the 
U/Ne spectrum. We then calibrated the U/ Ne spectrum using the recently-published 



uranium line list of iRedman et al.1 (|2011r) . which is derived from high-resolution 



Fourier transform spectrometer measurements. These two independent calibrations 
allowed us to easily identify emission lines in the hollow cathode lamp that do not 
correspond to known (classified) lines of either uranium or neon, and to compare 
the achievable precision of each source. Our frequency comb precision was limited 
by modal noise and detector effects, while the U/Ne precision was limited primarily 
by the signal-to-noise ratio (S/N) of the observed emission lines and our ability to 
model blended lines. The standard deviation in the dispersion solution residuals from 
the S/N-limited U/Ne hollow cathode lamp were 50% larger than the standard devia- 
tion of the dispersion solution residuals from the modal-noise-limited laser frequency 
comb. We advocate the use of U/Ne lamps for precision calibration of near-infrared 
spectrographs, and this H-band atlas makes these lamps significantly easier to use for 
wavelength calibration. 

Subject headings: atomic data — infrared: general — instrumentation: spectrographs 
— methods: laboratory — standards — techniques: radial velocities 
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1. Introduction 



The growing interest in high-precision astrophysical spectroscopy (such as the detection of 
low-mass exoplanets) requires the measurement of radial velocities at the meter per second (m/s) 
level. This corresponds to a Doppler shift in the astrophysical spectrum of 3 parts in 10 9 . Such 
precision is being achieved now in the optical with iodine cells or hollow-catho de lamps (HCLs) 
with s table heavy eleme nts such as thorium ( 232 Th) and a fill gas (usually argon) (IMarcy & Butler 
199ilMavoretal.l2009l). In the nea r-infrared (NIR), su b- 10 m/s precision calibration sources such 
as telluric C0 2 lines dFigueira et alJboioh and gas cells JlVIahadevan & Gell2009l : lBean et alJboid) 
are limited to narrow spectral regions, and thorium HCLs have few bright, calibrated lines com- 
pared to the visible. We propose that ur anium neon (U/Ne) HCLs are a more su itable calibration 
emission source in the NIR than Th/Ar ([Redman et al]|2011uRamsey et al.l l2008). 



The ideal calibration source would be a picket fence of bright emission lines or deep ab- 
sorption lines covering a wide wavelength coverage and exhibiting a uniform density and uniform 
intensities throughout the observed spectral region. Such a calibration source should not be too 
dense, because blended lines are difficult — if not impossible — to use for precise dispersion solu- 
tions. The calibration source should also be repeatable over long temporal baselines, and it should 
be relatively inexpensive and easy to operate. 

Laser-frequency combs (hereafter, LFCs) provide a picket fence of lines and are an ideal 
calibration source. Even with their high cost and technical complexity, they have quickly been 
adopted by the a strophysical communit y, which is anxious for the development of such a high- 



precision source. ISteinmetz et al.1 (|2008|) used the Vacuum Tower Telescope (VTT) solar telescope 



to demonstrate a comb with sunlight. Efforts are being made to develop and test such astro-combs 
by many groups in the optical for the High- Accuracy Radial velocity Planet Searcher (HARPS), the 
Trans-Atlantic Exoplanet Survey ( TrES), and the High- Accuracy Radial velocity Planet Searcher 
North (HARPS-N) spectrographs Li et alJbooa Iszentgvorgvi et alJl200a Iwilken et alJboid) as 
well as future European Southern Observatory (ESO) facilities like the Ech elle Spectrog raph for 
Rocky Exoplanet and Stable Spectroscopic Observations ( ESPRESSO) (IPepe et alJl2010D and the 
COsmic Dynamics and EXo-earth experiment (CODEX ) dPasauini et al. 2005 ). Our own efforts 



have focused on developing NIR laser-frequency combs (Ouinlan et al 



to enable the detec- 



tion of terrestrial-mass planets in the habitable zones of M dwarfs (the smallest and most populous 
stars in the galaxy) when coupled to the next generation of NIR precision radial velocity spectro 



graph s like the Habitable Zone Planet Finder (HZPF) that are now being built (IMahadevan et al. 
2O10h . 



In the meantime, HCLs remain the next best solution for the high-precision calibration of 
fiber-fed spectrographs in the NIR. Thorium is an ideal atomic emission source, as it has only one 
natural isotope, zero nuclear spin, a long half-life, and an abundance of lines. Uranium shares 
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all of these traits except the first: there are three naturally-occurring isotopes of uranium, but the 
abundance is dominated by 238 U, which makes up 99.275 % of the element, followed by 0.720 % 
of 235 Uand 0.005 % of 234 U0. Uranium also has approximately twice as many emission lines in 
the NIR (intrinsical ly). Since only a fraction of the thorium lines in the NIR have been published 
(|Kerber et al.ll2008l) . a more co mplete uranium lin e list represents a significant gain in the number 
of measured calibration lines. Palme r et al. dl98Ch published an atlas of uranium for the optical, 
and Englemanl (J2003) commented on the significantly higher density of bright lines in the emission 



spectr a of uranium compared to thorium. There ar e now three uranium line lists: [Palmer et al. 



(I1980|) . which covers the optical up through 910 nm. lRedman et al] (1201 II) (he reafter. Rl 1), which 
features uranium lines between 850 nm and 4000 nm, and Conway et al] ( 1984 ). whi ch covers the 



NIR between 1800 nm and 5500 nm. We refer the interested reader to the work of Kerber et al 



(I2007|) for a review of the development of these lamps, and to Kerber et al] (I2008|) for a detailed 
NIR thorium line list. 

Bright lines from the fill gases of hollow-cathode lamps are useful for moderate precision 
wavelength calibration, but are a hindrance for sensitive high -resolution instrumen ts that require 
a very-high-precision calibration and minimal scattered light. IWhaling et al] (|2002|) have demon- 
strated that pressure changes shift the measured wavenumbers of transitions from the 4f, 5g, and 
6g levels of Ar I. As such, these lines are not suitable for precision radial velocity measurements 
below tens of m/s (|Lovis & Pepdl2007|) . In the NIR, bright argon lines hamper the analysis of the 
fainter Actinide lines that are more stable and suitable for wavelength calibration. Neon, on the 
other hand, has relatively few lines beyond 900 nm, and thus produces less scattered light and pixel 
saturation — see figure 1 of Rl 1 for a comparison between Th/Ar, U/Ar, and U/Ne. These factors 
make U/Ne a potentially useful calibration source for precise NIR spectroscopy. 

Spectral atlases are essential to the initial wavelength calibration of grating spectra. They 
provide both approximate relative intensity and wavelength information that the human brain can 
quickly compare to an observed spectrum in a fraction of the time that would otherwise require ex- 
tensive computer programming and processing time. Atlases are useful both for gross wavelength 
calibration of an unknown spectral region (e.g., when a spectrograph is commissioned), and for 
the identification of specific spectral regions (e.g., when a telescope user obtains a spectrum at a 
particular grating setting). Once the observed spectral region has been identified, the user can refer 
to the original line list for specific line information, and measure the line centroids to extract the 
dispersion solution. 

We have created an atlas of uranium-neon between 1454 nm and 1638 nm using a s imul- 
taneous laser-frequency comb (|Murphy et al]|2007l : lOsterman et al]|2007l : iBraje et al] 120081) with 



^http : / /physics . nist . gov/ PhysRef Pat a /Handbook /Table s/uraniumt able 1 . htm| 
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a second calibration fiber to provide a first-order estimate of the wavelength calibration of the 
hollow-cathode spectrum. We then refined the dispersion solution using the uranium line list of 
Rl 1. We observed this spectrum at a current (14 mA) and resolution (50 000) that are more typical 
of astrophysical spectrographs than the h igh-curren t (75 mA to 300 mA), high-resolution Kitt Peak 
Fourier transform spectrometer (FTS) ( Braultfl976 ) spectra that were used to measure the uranium 
wavelengths in Rl 1. The format of this atlas makes it suitable for day-to-day use by astronomers 
both for visual inspection and for wavelength calibration, and its layout is comparable to Th/Ar 
atlases that are in routine use at telescopes around the worlc§. 

A number of stable fiber-fed instruments that are now being planned or are in the design 
stages could benefit from the atlas published here. The Sloan Digital Sky Survey III (SDSS-III) 
Apache Point Observatory Galactice Evolution Experiment (APOGEE) galactic kinematics survey 
is currently commissioning a 300-fiber R^ 22 5 00 H-band fiber-fe d spectrograph covering the 
spectral region between 1530 nm and 1680 nm (IWilson et alj|2010|) . Ne w high-resolution NIR 
spectrographs like the Immersion Grating Echelle Spectrograph (iSHELL) (ITokunaga et al]|2008|) 
(already under construction) can benefit from better calibration sources. The next generation of 
high-resolution NIR spectrog raphs built specifically for precision RV surveys, like the Habitable 
Zone Planet Finder (HZPF) (IMahadevan et al.ll2010|) . the Calar Alto high-Resolution search for 
M dwarfs with Exoearth s with Near-infrared and optical Echelle S pectr ographs (C ARMENES ) 
( Ouirrenbach et al. I2OI0L and SpectroPolarimetre InfraROUge (SPIRou) ( Donati in preparation ), 
are now in their planning stage and the U/Ne calibration source we present here may help them 
achieve their desired precision. 

We describe the spectrometer and our H-band frequency comb in §2. In §3 we briefly describe 
how we obtained our spectra. In §4 we explain how we reduced and wavelength-calibrated the LFC 
and U/Ne data. In §5 we discuss some of the limitations of the atlas and compare it to the precision 
of the frequency comb. This atlas is available online at URL from publisher. 



2. Instrumentation & Calibration 



2.1. The Fiber-Fed NIR Pathfinder Spectrograph 



The Pathfinder testbed used for this work is similar to the instrument we used as a testbed 



to demonstrate 7 m/s to 10 m/s NIR RV precision on sunlight in Ram sey et al.l (|2008). Details 
of our current setup, and improvements, are described in iRamsey et al. J2010h . but we mention 
key design elements here for completeness. Pathfinder is a fiber-fed warm-pupil cross-dispersed 



, noao . edu/kpno/specatlas/thar/thar . html| and |het .as . utexas . edu / HET /net web / Instrument s/K 
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echelle spectrograph that yields a spectral resolution A/ A A ~ 50 000. The input fiber coupling 
is designed to match the f/4.2 output of the fiber from the Hobby-Eberly Telescope (HET) fiber 
instrument feed. A pair of commercial NIR achromats optimized for the 0.75 /mi to 1 .55 ftm region 
image the beam from the HET and calibration fibers onto a 100-/xm slit. The diverging beam is 
folded by a small gold mirror on the way to the collimator, which in turn directs the light to a 7 1 .5° 
blaze angle (R3) echelle grating. While this grating is undersized in the dispers i on dir ection for 
the 100-mm Pathfinder beam, it is twice as efficient as that used in lRamsey et al.l (|2008l) . This R3 
echelle operates in-plane at an off-Littrow angle of « 5° to achieve high dispersion. This mode 
of operation (no 7 angle tilt of the echelle) ensures that there are no slit curvature effects. A 
gold-coated 150-1/m m grating blazed at 5.46° provides cross dispersion, enabling our 1024 x 1024 
HAWAII- IK array (|Hodapp et al.ll 19961 12004|) to simultaneously image 7 to 8 echelle orders in the 
Y band and 3 to 4 ord ers in the H band. We use the same camera system as described in detail 



111 



Ramsey et al.1 (|2008|) . which uses a parabolic mirror with a weak coma-correcting lens near the 



entrance of the NIR dewar. The uncooled Pathfinder testbed is made possible by liquid-nitrogen- 
cooled (78 K) thermal blocking filters that suppress the out-of-band thermal radiation to a high 
degree (fMahadevan et al.llin preparation!) . 



In its H-band configuration, Pathfinder has three optical fibers (one stellar fiber leading to 
the focal plane of the HET and two calibration fibers that lead to our calibration bench) that can 
simultaneously be illuminated. These fibers were polished and assembled on-site at the HET and 
then aligned with a 100-/mi slit, shown in Figured] The alignment of the three fibers with the slit is 
not perfect and this causes a constant offset in the wavelength calibration from the three fibers. By 
feeding laser-frequency comb light through all three fibers simultaneously, the dispersion solution 
of each fiber can be independently measured to obtain the offset between the different fibers (see 

§EB. 



2.2. The H-band Frequency Comb as a NIR Calibration Source 

The National Institute of Standards and Technology (NIST) and University of Colorado at 
Boulder (CU) H-band laser- frequency comb is built around a 250 MHz passively mode-locked 
Erbium fiber laser. The laser is frequency-stabilized by locking both the repetition rat e f rev (direct 



detec tion) and the carrier-envelope offset frequency / ceo if - 2f optical heterodyne, IJones et al. 



2000|) to a Rb clock that is steered on the long term by the global positioning system (GPS). In this 



configuration, the passively mode-locked laser generates modes at frequencies of 



f n — n x f TC p + f ce0l 



(1) 
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centered at 1 550 nm with a bandwidth of about 70 nm. The accuracy of the mode frequen- 
cies is limited by th e accuracy of the GPS -stabilized Rb clock, and is better than 1 part in 10 11 
(|Ouinlan et al.ll2010h . For use with Pathfinder, the mode spacing is increased using Fabry-Perot 
filter cavities. Two identical cavities are employed, with 25 GHz free spectral range and finesse of 
about 2000; one is placed immediately after the mode-locked laser, and the other is placed after 
a series of optical amplifiers, see Figure [21 The lengths of the cavities are actively locked by use 
of piezo-electric actuators to the laser-frequency comb for maximum transmission of one set of 25 
GHz-spaced comb modes. After filtering and amplification, the 25 GHz comb is broadened in a 
highly nonlinear fiber, then transmitted via a single-mode fiber to an integrating sphere, where a 
fraction of this light enters a 300-/im fiber that leads to the spectrograph. After filtering, the comb 
equation for the 25 GHz comb is 



/„ = (n + 100n)x/ rep +/ c 



(2) 



where n Q is the index of one comb mode transmitted through the filter cavity, measured with a 
wavemeter. For these experiments, n was 782971, / rep was 250.1132048 MHz, and / ceo was 
70.11320482 MHz. Fur ther details of the comb and our tests with Pathfinder are presented in 
Yeas et al.1 din preparation!) . 



3. Experimental Design & Observations 

The simultaneous frequency comb and U/Ne spectrum for this atlas were collected as part of 
an on- sky commissioning run of the Pathfinder spectrograph at the HET. Light from the frequency 
comb was fed into the primary calibration fiber, while U/Ne light was fed simultaneously through 
the secondary calibration fiber (the HET stellar fiber was not needed during this experiment). We 
ran a Photron§ U/Ne hollow-cathode lamp (Part number P863, Serial Number HGL0170) continu- 
ously at 14 mA. The Pathfinder HAWAII IK detector can only cover about 20% of the observable 
spectral range, so the grating settings had to be adjusted incrementally throughout the experiment. 
At each grating setting, we obtained ten 30-second U/Ne and frequency-comb images and five 30- 
second flat-field images from an quartz tungsten halogen bulb. All exposures were the same length 
so that we could take exposures continuously. 30-second dark frames were taken throughout the 
observing run and used to background- subtract the exposures. The flat-field exposures were taken 
with neutral density filters in front of each fiber, which had to be removed during the U/Ne and 



Identification of specific products is for scientific information only and does not constitute an endorsement by 
NIST. 
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comb exposures. During these exchanges, a shutter in front of the fiber output in the pathfinder 
experiment was closed to minimize the effects of persistence in our images. The frequency comb 
was monitored continuously to ensure the stability of n , f rep , and f ceo . 

An annotated map of U/Ne atlas, created by overlapping adjacent spectral images, can be 
seen in Figure [31 In each order, the frequency comb spectrum is above the corresponding U/Ne 
spectrum. Wavelength increases to the left and down. The index of every tenth frequency comb line 
has been marked, and the frequency of each comb line can be found by using Equation^ which is 
also printed on the image for easy reference. The blackbody radiation of the hollow-cathode lamp 
is visible as a faint continuum throughout the U/Ne spectrum; this continuum is not present in the 
frequency comb spectrum since the comb is not a thermal source. Several bright neon lines have 
been marked throughout the map. 



4. Analysis 
4.1. Data Reduction 

We constructed a median-averaged flat-field frame and a median-averaged U/Ne and laser- 
frequency comb image for each grating setting (in order to eliminate cosmic rays), as well as a 
median-averaged image of our darks to form a master dark. This master dark frame was sub- 
tracted from each of the flat field and science frames. We used subroutines from REDUCE 



(|Piskunov & Vale nti 2002) to model and normalize our master flats, subtract scattered light from 
the U/Ne/comb images, and extract the U/Ne and comb spectra. After wavelength-calibrating the 
spectra (see § 14.21) , the U/Ne and comb spectra at each grating setting were resampled and averaged 
together to produce a continuous spectrum from 1454 nm to 1638 nm. Outside of this range, the 
comb lines were too faint to identify. 



4.2. Wavelength Calibration 

4.2.1. Wavelength Calibration of the Frequency Comb Spectrum 



Wavelength identification of the frequency comb spectrum was conducted by manually iden- 
tifying the comb line indices (n) spe cific to a given g rating setting and fitting the comb peaks with 
Gaussians using MPFIT routines (IMarkwardtl l2009h in Interactive Data Language (IDL). Unlike 
a hollow-cathode emission spectrum, frequency comb lines are evenly spaced and identical in ap- 
pearance. Because n , f re p, and / ceo are known, we can determine n by feeding a single fiber with 
both comb light and a calibration with a well-known spectrum. In our case, since we know the 
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fibers are separated by less than a pixel in dispersion 1 GHz), while the comb lines are sepa- 
rated by about 25 pixels 25 GHz), we used bright (but not saturated) uranium and neon lines 
from the adjacent fiber to determine the index of the nearest comb lines and thereby determining 
the index (and frequency) for all lines in that order. Note that the determination of n requires much 
less precision (by two orders of magnitude) than our later comparison between the FTS uranium 
wavelengths and the comb-based uranium wavelengths. 

The LFC dispersion solutions were found by fitting the comb line wavelengths as a function 
of modeled Gaussian centroids with a fourth-order polynomial, ignoring outliers that fell more 
than three standard deviations from the mean residual to avoid the influence of coincident cosmic 
rays or stray light from bright lines in adjacent orders. An example comb spectrum can be seen 
in panel A of Figure HI The residuals of this dispersion solution can be seen in the same figure as 
diamonds; the scale for the residuals are on the right-hand axis. 

In the absence of systematic sources of noise, the RMS of the frequency comb dispersion 
solution residuals should decrease as the square root of the signal-to-noise-ratio (hereafter, S/N) 
of the lines that were used to solve the dispersion solution. However, we found that the frequency 
comb RMS residuals plateaued at about 0.000 4 nm. We hypothesize that this noise floor is induced 
by a combination of modal noise and detector inhomogeneities. 



Fiber modal noise ( Baudrand & Walked 1200 II) is an inevitable result of the finite number of 



modes propagating in the fiber, and leads to additional noise in all fiber-fed spectrographs. The 
modal noise in NIR fiber-fed instruments is exacerbated (compared to the visible) because the 
longer NIR wavelengths lead to fewer modes in the optical fiber. In spite of mitigation efforts, 
like putting the LFC light through the integrating sphere and fiber agitation, we discovered that 
the primary and secondary calibration fibers were still impacted by modal noise. While the effects 
of this were mitigated for on-sky observations by averaging together many cumulative exposures, 
the modal noise does affect the U/Ne and LFC calibration experiment at the level of a few tens of 
m/s, and sets a floor on the achievable precision. Better agitation and scrambling mechanisms are 
expected to mitigate this problem. 

The uncertainties of the frequency comb lines are the sum in quadrature of two components: 
the statistical uncertainty of the line (as measured by the width Wf com bj and S/N S/Nf Com bj of the 
model Gaussian), and the uncertainty due to modal noise and detector inhomogeneities: 




2 

+ (0.000 4 nm) 2 (3) 



Technically, we should also include the uncertainty in the frequency comb wavelength, but this is 
several orders of magnitude smaller than the other uncertainties and therefore negligible for our 
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purposes. 



4.2.2. Wavelength Calibration of the Uranium-Neon Spectrum 

The measurement of frequency comb centroids and the calculation of the dispersion solution 
from their corresponding frequencies is a simple task compared to the wavelength calibration of a 
hollow cathode spectrum because the latter spectrum contains blended lines. The dispersion solu- 
tion of uranium neon was first estimated by offsetting the frequency-comb dispersion solution of 
the same grating order by 0.005 nm, the approximate separation between the fibers. This separa- 
tion varies across the order slightly, and is different for each order and for each grating setting, so 
the dispersion solutions had to be solved using the U/Ne lines alone. In an effort to derive the most 
accurate and precise dispersion solutions possible, we solved our U/Ne dispersion solutions in a 
consistent and automated manner. The following process is applicable to the measurement of any 
hollow-cathode dispersion solution. 

Each U/Ne spectrum was modeled as the sum of Gaussians, one for each line in the uranium 



and neon line lists of Rl 1 and lSaloman & Sansonett i (2004). For the uranium lines, we used the 



relative intensities of Rl 1 to guide our initial guesses. We also added a constant offset to the model 
for the c ontinuum. This m odel spectrum was iteratively modified to fit the observed spectrum using 
MPFIT iMarkwardl2009L We then restricted the dispersion solution to a careful selection of lines. 



First, we neglected to use any neon line s for our dispersion solutions. Some of the gas lines of 



argon are susceptible to pressure shifts ([Whaling et al.ll2002 ): while it is not known whether or not 



any of the neon energy levels are susceptible to these shifts, we thought it prudent to ignore these 
lines from the beginning. Second, we restricted ourselves to those lines with a S/N above 100. We 
lowered this threshold if fewer than 10 lines per order had a S/N above 100, but the threshold was 
usually above 40. Third, we used the modeled Gaussians goodness-of-fit to guide our selection. 
Lines that did not appear in our observed spectra had modeled widths of zero, and lines that were 
poorly modeled exhibited relatively large or noisy residuals. In an effort to avoid the influence of 
heavily -blended and poorly-modeled lines, we calculated the residuals as the difference between 
the observed spectrum and the modeled line (neglecting the other modeled lines). If either the 
mean difference or standard deviation in the difference of the three closest pixels was larger than 
5% of the line peak, the line was not used to solve the dispersion solution. This threshold was 
increased in the event that there were fewer than 10 usable lines, but was usually below 8%. The 
uranium lines that met these criteria accounted for approximately 20% of the uranium lines from 
Rl 1 in this spectral range, and are marked throughout the atlas with an asterisk ('*'). 



An example uranium spectrum and model (in grey; in red in electronic version) is shown in 
panel B of Figure @J Using the above-mentioned selection of uranium lines, we solved the dis- 
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persion solutions by fitting a fourth-order polynomial to the uranium wavelengths as a function 
of the modeled line centroids. The residuals to the dispersion solution are shown as diamonds. 
The uncertainties of the dispersion solution residuals are the sum in quadrature of three compo- 
nents: the statistical uncertainty of the model line, the uncertainty due to modal noise and detector 
inhomogeneities, and the uncertainty of the wavelength of the line, as measured in Rl 1: 



v W(fc) 2+<aMo4nm)2+ *'< (4) 

where W v / Ne j is the width of the model line and S/Nu/mj 1S tne signal-to-noise ratio of the model 
line. The residuals from the fit are shown in panel C of Figure @J 



5. Results 

The extracted uranium-neon H-band atlas is available online; an example of this atlas is pre- 
sented in Figu re [51 The wavelengths and u ncertainties of the uranium and neon lines are taken 



from Rll and iSaloman & Sansonettil (120041) . respectively. Lines of an unknown origin (marked 
with an '?') are also from Rl 1, and are most likely uranium lines with unknown upper energy lev- 
els. The flux has been normalized to the brightest feature in the spectrum, the (saturated) neon line 
at 1523.487 65 nm. An electronic version of this atlas, and the wavelength-calibrated spectra used 
to create it are available at (URL for VIZeR, from publisher). This atlas provides astronomers 
with a useful visual reference for deciding upon which spectral region to use for calibration, as 
well as a means of confirming the observed spectral region. 

There are many features of note throughout the atlas, and we highlight a few of them here: 

1. There is a uranium line at 1511.052 33 nm that falls within 0.006 nm of a known neon 
line. The uranium line is based upon its existence in U/Ar spectra. Based upon the relative 
photon fluxes from Rl 1, and the fluxes of the nearby uranium lines, both lines contribute an 
approximately equal amount of flux to the feature. 

2. Some of the many blended lines are readily apparent throughout the spectrum, and confirmed 
both by the FTS spectra and the Pathfinder spectrum (e.g., the pair of uranium lines at 1634. 1 
nm). 

3. There are several unidentified lines throughout the spectrum, such as the line near 1555. 1 
nm. These lines are not pres ent in the uranium line list of Rl 1, nor in the neon line list of 



Saloman & Sansonettil (|2004). As such, they are not marked in the atlas, but are apparent 
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with a visual inspection. Others include the lines at 1578.6 nm and 1480.6 nm. These 
might be unknown uranium or neon lines, but they may also be due to contamination in our 
lamp. No known lines of any ele ment were found at these locations in the Atomic Spectral 
Database dRalchenko et alJboi II) . 



4. Sharp artificial features sometimes appear in the spectrum when bright lines near the edge 
of the detector shine scattered light onto nearby grating settings. These illumination patterns 
produce discontinuities in the spectral averaging algorithm; an example can be seen near the 
neon line at 1493.388 63 nm. 



5.1. Comparison Between the Wavelength Calibrations 

A comparison between the precision of the LFC and U/Ne hollow cathode dispersion solutions 
is shown in the top two panels of Figure [61 The standard deviation of the LFC residuals was 
0.000 8 nm, while the standard deviation of the U/Ne hollow cathode residuals was 0.001 2 nm. 
The former was limited by modal noise and detector imperfections, while the latter was limited by 
blended lines in the spectrum, which offset the modeled centroids in a systematic fashion. 

The limitation induced by modal noise was largely unexpected, as we were agitating the 
output end of the fiber, whi ch is known to reduce modal noise in the visible and very near-infrared 



(|Baudrand & Walkerll2001|) . We suspect that this modal noise was caused primarily by the way 
in which we fed light into the fibers: the optical alignment of the calibration system was done by 
eye (e.g., in the optical), and it's possible that different indices of refraction from the beam splitter 
caused a slight offset in the position of the H-band light relative to the optical. By injecting fewer 
modes into the fiber, we diminished the number of modes at the output. We anticipate that we can 
reduce modal noise in the H-band in future Pathfinder experiments by aligning our fibers in the 
NIR and using a double scrambler. 

The uranium hollow-cathode lamp is dense with emission lines, especially in the NIR. This 
density is a double-edged sword, as it provides many potential calibration lines throughout the 
NIR, but also a large number of blends, many of which cannot be used for precise radial velocity 
measurements because they obfuscate the true dispersion solution. The precision of a given line 
depends upon the precise nature of the blend — two equally bright lines will induce more of a 
shift in their measured line centers than will a blend where one line is significantly weaker than 
the other. This is a significant advantage that a frequency comb has over hollow-cathode emission 
lamp sources, since blends in the former are completely avoidable, if the comb lines are optimally 
filtered. 



Unlike measuring the precision of a wavelength calibration source, measuring its accuracy 
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requires a comparison to wavelengths to an independent standard. This is particularly important 
in the case of a hollow cathode spectrum, where the observer does not always know if there are 
unknown, weak lines introducing systematic errors into the dispersion solution. This is undoubt- 
edly the reason that the thorium-argon lamps on HARPS failed to reveal the underlying detecto r 
stitching that was so readily seen when the detector was illuminated by a LFC (Wi lken et al. 



In an effort to estimate the accuracy of our U/Ne dispersion solutions, we measured the wave- 
length difference between the calibration fibers by feeding them simultaneously with frequency 
comb light. The dispersion solution differences between these fibers for each of three different 
orders is shown in Figure [7J (solid curves). Also plotted on this figure is the difference between 
the frequency comb and the U/Ne dispersion solution at the same grating setting (dashed curves). 
In almost all cases, the U/Ne dispersion solution matches the frequency comb dispersion solution 
to within the uncertainties of the uranium lines. Unfortunately, we could not perform this analy- 
sis at every grating angle for lack of time, but the bottom panel of Figure [6] shows the difference 
between U/Ne and frequency comb dispersion solutions throughout the H-band. The differences 
in the three orders of Figure [7J are highlighted in grey (green in electronic version). Across the 
H-band, the dispersion solutions exhibit approximately the expected fiber separation. 

It is possible that one could achieve better precisions with a U/Ne hollow cathode source 
through a more careful selection of emission lines, but this experiment provides a unique oppor- 
tunity to assess both the accuracy and precision of the chosen lines, which would not be possible 
without the LFC (or another highly-accurate and highly-precise fiducial). Experiments at lower 
precision should use a more conservative selection of lines, since blending will be more frequent, 
while experiments at a higher precision should be able to use a larger selection of uranium lines. 



6. Conclusions 



Uranium hollow-cathode lamps have a significantly larger number of usable bright lines than 
thorium hollow-cathode lamps in the NIR Y, J and H bands. This high density of lines enables pre- 
cise wavelength calibration for the vast majority of NIR spectroscopic applications. This atlas of 
U/Ne significantly reduces the barriers to using such lamps in modern NIR spectrographs by pro- 
viding an easy-to-use visual reference. In addition, the identification of the least-blended uranium 
lines enables the user to avoid regions and lines that are unsuitable. This work also serves to il- 
lustrate the immense promise and advantage of LFCs for precision astrophysics. While frequency 
combs offer the best solution for ultra-precise wavelength calibration, their cost and complex- 
ity currently limits their applications to astronomical instruments that require an extremely high 
level of wavelength calibration precision over multi-year baselines. Relatively inexpensive hollow- 
cathode lamps are the most cost-effective approach for the majority of astronomical spectrographs. 
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Fig. 1. — Top: Pathfinder Spectrograph input fiber feed. The three 300-/im fibers, separated by 
125-/im (outer diameter) fibers that acted as spacers for the larger fibers. This fiber pseudo-slit was 
assembled and polished on-site. From left to right, these are the secondary calibration, primary 
calibration, and science (stellar) fibers. Bottom: The same fibers, re-imaged onto a 100-/im slit. 
The slit is shown in the foreground, covering the top and bottom of the fibers. In this image, the 
primary calibration fiber is transmitting light. The slit was aligned to optimize the orientation of 
the primary calibration and HET fibers. Note that the HET fiber was not used for the experiment 
described in this paper. 
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Fig. 2. — Schematic of 25 GHz frequency comb. A 250 MHz mode-locked Er: Fiber laser is fre- 
quency controlled by locking both repetition (f rep ) and carrier-envelope offset (f CEO )frequencies. 
Light from the 250 MHz frequency comb is coupled into a 25 GHz FSR Fabry-Perot cavity com- 
prised of Mi and M 2 with lenses Li-L 4 matching the modes of the transmission fiber to the cavities. 
A piezoelectric actuator (PZT) controlled by a servo is used to actively lock the cavity length to the 
frequency comb. The 25 GHz comb is optically amplified, filtered again in an identical Fabry-Perot 
cavity, and then spectrally broadened in highly nonlinear fiber (HNLF.) Standard single-mode fiber 
delivers the broadened 25 GHz comb to the spectrograph, where it is launched into an integrating 
sphere. 



- 19- 




It I I ■>• XI 'It > 1 1 * I 



OrdorJO f 1499.04t 60 nm 

-60 .50 -40 -30 -20 

I , J 1 I i 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



1 1 . 1 


1 1 1 1 1 > II 


1 1 1 1 1 1 1 > > 1 


■ 1 1 1 1 




-10 o 

1 1 


10 20 
1 1 


30 

m 




40 



I I I I I II I I I I M I I I I I I I I 1 I I 



Oi*r39 f 1541.180 32 nm 

-250 -250 -240 



1 1 



1 



■220 -210 " ?QD -190 -160 -170 

i I I I 

lil M I H 1 1 1 




Orter38 
-(SO 

l I I 



1 1 1 I i J J i I 

I I I! I l l I I I M M i n I I I I I M M 1 1 I HUH MM I » M > Ml M I I M | M '| M 1 M 1 U I J 1 1 1 I I M «■«> ■■ j J • 




•1523.467 65 rim 

-160 150 -140 



-450 



-350 -340 
i I 



440 -430 -420 -410 -400 -390 "380 370 -360 

I 1 1 i I I 1 i 

, null IIIIMIMIIHIIMI mmIMMIMMIIIMMHIMMIIMIIIIM Ml HimilimiMMIHI 

t 1602.71463 nm 

Frequency comb calculation (201 August 1 2] 

(]rsl - (7B2971 + 100'rs)'250.1 132040 + 70.1 132Q4S2 MHz " 550 "540 

1 I 



* 4 i i i 

Order 37 




SO 290 300 31 320 330 340 350 360 370 



i i i t r i ; * { i 

t»..»l . « , » .( It* »t i i M «i i M<* f. MI 1*41 I «* ' • •**»• MW' 

90 100 




110 



120 



130 140 



II* II 

190 



150 160 170 18° 

l I £ 1 i 1 1 1 1 I J 1 f 

[,,,,, iiiillll Ml Mil Mill I I I M M M I M >l I I I l I I . I l l I I I > I I I H I I M I I M 1 t i • i ■MMMIMI 'j"M 'll 

IM^H .*»►•«». ■ «t >«4"li<*"**'Mi '• • * H tw^mHWI 

t 1 499.041 50 nm 

-120 -110 -100 -90 -80 -70 -60 -50 -40 '3° " 10 

1 i J 1 i I 1 1 i I 

Hi I I . I I i I I I II I I I I I II I I I I I I I I M I I I I I * I I I I I I I l I I < I > I I I II I I II I I 

- *4»<* * 



Mil ii Htii i I I' 

-330 -320 -310 " 300 

I I i i 

1 1 1 1 1 1 1 M I H 1 1 M 1 1 1 1 1 1 1 1 II I M (1 1 M ' i ■ ' 

I i i. 4 ■ i « I' I ■ i 



1 1 M I I I I 1 1 I 1 1 1 1 1 M 11 1 1 H 1 1 I I • > 
I 1 1 i tut »• II it • ♦ i I * 



t 1541.180 32 nm 



■290 



•280 
1 



-270 
J 



-260 
I 



-530 
I 



-520 
i 



-510 



' 1 

-5D0 
1 



-490 
I 



I l i|t 

-480 
i 



, ... i | ll 



i i f { 
huh n i n 

». < .|» I t • 



-470 
1 



-460 



-450 -440 



-430 
1 



, | IIIMMII IIIIIMM I IMIIHIIHMin Mill,. 

,,,,,, « h i » I • i 1 ' 1 I ■ • " * ' • » 1 * * " 

t 1640.973 72 nm 



Fig. 3. — A map of uranium-neon, taken with the Pathfinder spectrograph, and presented in two 
parts. The orders are indicated in the top half of the figure. The integers indicate the frequency 
comb index of the indicated line, which can be used with the frequency comb equation to cal- 
culate the frequency of any given comb line: /„ = (n +l00n) x f rep + f ceo = (782971 + lOOn) x 
250.1132048 MHz + 70. 1 1320482 MHz. Five bright n eon lines have been indicated throughout 
using the neon line list of ISaloman & Sansonettil (|2004l) . 
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Fig. 4. — An example comparison between the laser-frequency comb dispersion solution residuals 
and the uranium- neon dispersion solution residuals. A: The centroid of each frequency comb line 
was found by fitting the spectral features with independent Gaussians (not shown). The wavelength 
as a function of Gaussian centroid was fit with a fourth-order polynomial. The residuals of this 
dispersion solution fit are over-plotted as diamonds on the frequency comb spectrum. The scale 
for the residuals are on the right-hand axis. B: The observed U/Ne spectrum was m odeled as the 
sum of Gaus sians, one for each line in th e uranium line list of iRedman et al.1 (|2011r) and the neon 
line list of ISaloman & Sansonettil (120041) . The best-fit model, using a chi-squared minimization 
routine, is shown here in grey (red in electronic version). Only certain model lines were used to 
solve the dispersion solutions (see § |4.2.2| for details). The dispersion solution was fit with a fourth- 
order polynomial, and the residuals are over-plotted as diamonds. C: The difference between the 
observed U/Ne spectrum and the model U/Ne spectrum. Note that we saw several emission lines 
that do not correspond to known uranium or neon lines in this spectrum. These lines (e.g., the line 
near 1621.7 nm, and at least two lines near 1623.5 nm) were not well-modeled, and show large 
residuals. 
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Fig. 5. — An example image from the uranium-neon atlas. The entire atlas is available online. 
Each spectrum has been wavelength-calibrated by the brightest uranium li nes. The wave l engths 
and uncertainties of the uranium, unknown, and neon lines are taken from iRedman et al~ ( 2011 ) 



and lSaloman & Sansonettil (|2004). Those lines used for the uranium-modeled dispersion solutions 
are marked with an asterisk ('*'). The flux has been normalized to the brightest feature in the 
spectrum, the (saturated) neon line at 1523.487 65 nm. 
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Fig. 6. — A comparison between the precision of the frequency comb and U/Ne dispersion so- 
lutions, and an estimate of the accuracy of the U/Ne dispersion solutions. Top: The dispersion 
solution residuals from fitting each comb line with a Gaussian. The standard deviation of these 
points is 0.000 8 nm. Middle: The dispersion solution residuals from the U/Ne spectrum, where 
each FTS-measured uranium and neon line has been modeled as a Gaussian and fit to the observed 
spectrum. Only the brightest, best-fit uranium lines of each order were used to calculate the dis- 
persion solutions. This method allows us to use blended lines to a much higher degree of precision 
than the simpler peak-fitting algorithm. The lines used in the bottom plot are marked throughout 
the atlas with an asterisk ('*'). The standard deviation of these points is 0.001 2 nm. The uncer- 
tainties are larger beyond about 1610 nm because the uranium lines in this region are weaker (in 
part because of a lower spectral response). Bottom: The difference between the U/Ne and comb 
dispersion solutions across the H-band. The differences presented in Figure [7] are shown here in 
grey (green in electronic version). The mean offset between the fibers is 0.005 1 nm, and the single 
standard deviation of these offsets is 0.001 3 nm. 
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Fig. 7. — The fiber separation between the primary calibration and secondary calibration fibers for 
three different echelle orders. The solid line indicates the mean dispersion solution difference over 
52 images as measured when both fibers were simultaneously illuminated by laser frequency comb 
light. The uncertainties are the single standard deviation in the dispersion solution differences over 
these 52 images. The dashed line indicates the difference between the U/Ne spectrum and the 
laser frequency comb over a single image. The uncertainties on these points are the uncertainties 
in the uranium lines used to define the U/Ne dispersion solution. In general, the U/Ne dispersion 
solutions match the dispersion solutions of the comb, within the uncertainties of the uranium lines. 
Note how much the dispersion solution changes within a given order — this is the reason that 
dispersion solutions from one fiber cannot be applied directly to the corresponding order of a 
different fiber by a constant offset, without the loss of some accuracy. 



